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ABSTRACT 

Aims. We investigated the chemistry of nitrogen-containing species, principally isotopomers of CN, HCN, and HNC, in 
a sample of pre-protostellar cores. 

Methods. We used the IRAM 30 m telescope to measure the emission in rotational and hyperfine transitions of CN, 
HCN, ^^CN, H^^CN, HN^^C, and HC^'^N in L 1544, L 183, Oph D, L 1517B, L 310. The observations were made along 
axial cuts through the dust emission peak, at a number of regularly-spaced offset positions. The observations were 
reduced and analyzed to obtain the column densities, using the measurements of the less abundant isotopic variants 
in order to minimize the consequences of finite optical depths in the lines. The observations were compared with the 
predictions of a free-fall gravitational collapse model, which incorporates a non-equilibrium treatment of the relevant 
chemistry. 

Results. We found that CN, HCN, and HNC remain present in the gas phase at densities well above that at which CO 
depletes on to grains. The CNiHCN and the HNC:HCN abundance ratios are larger than unity in all the objects of our 
sample. Furthermore, there is no observational evidence for large variations of these ratios with increasing offset from 
the dust emission peak and hence with density. Whilst the differential freeze-out of CN and CO can be understood in 
terms of the current chemistry, the behaviour of the CN:HCN ratio is more difficult to explain. Models suggest that 
most nitrogen is not in the gas phase but may be locked in ices. Unambiguous conclusions require measurements of the 
rate coefficients of the key neutral-neutral reactions at low temperatures. 
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1. Introduction 

Observations of rotational transitions of molecules and rad- 
icals play a key role in deriving information on solar-mass 
objects in the early stages of gravitational collapse. The 
variation of the intensities of the emission lines can be in- 
terpreted in terms of the chemical reactions and gas-grain 
interactions occurring in the medium, and the line profiles 
and frequency shifts in terms of the kinematics of the col- 
lapsing gaseous material. Indeed, apart from infrared obser- 
vations of dust continuum emission, which yield no chemical 
or kinematical information, measurements of radio transi- 
tions of molecules provide the only means of probing the 
evolution of pre-protostellar cores. 

An obstacle to the use of molecular line emission to 
study the early stages of star formation is the propensity of 
some molecules to freeze on to the grains at the low temper- 
atures, T « 10 K, which prevail. Observations of prestel- 
lar cores have shown that the fractional abundances of the 
carbon-containing species, CO and CS, decrease strongly 
towards the core centres, whereas the fractions of the 
nitrogen-containing species, N2H"'" and NH3, either remain 
constant or even increase towards the centre, where the 
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density is highest (jTafalla et al.l[200^ . Differential freeze- 
out of the C- and N-containing species on to the grains 
was the generally-accepted explanation of these observa- 
tional results. However, recent observations of the NO rad- 
ical have demonstrated that the real situation is more com- 
plicated. A comparison of the profiles of NO and N2H'^ 
along cuts through the prcstellar cores L 1544 and L 183 
(Akyilmaz ct al. 2007], hereafter A07) has shown that the 
fractional abundance of NO, unlike that of N2H+, decreases 
towards the centres of these cores (their centres being iden- 
tified with the peak of the dust emission). Thus, not all 
nitrogen-containing species remain in the gas phase at den- 
sities > 10^ cm~^, which prevail in the central regions. 
On the other hand, st ill more recent observations of CN 
(|Hilv-Blant et all l2008[) have shown that the emission of 
this radical follows closely the dust emission in both L 1544 
and L 183. In so far as these two objects are representative 
of their class, it appears that the adsorption process must 
somehow differentiate between nitrogen-bearing species. 

In the present work, we extended our observations of 
nitrogen-containing species to include isotopomers of HCN. 
In addition to L 1544 and L 183, we have studied Oph D, 
L 1517B, and L 310. The observations are described and 
analyzed in Section [5] Sections [3] and S] describe the obser- 
vations and present estimates of N-bearing species abun- 
dances. In Section [5j we consider the chemical processes, 
including freeze-out on to the grains, which determine 
the gas-phase abundances of key N- and also C- and O- 
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Table 1. The sample of cores observed. Note that pOph D is also known as L 1696A. 



Source 


Q2000 


^2000 


Aa, Ar 


kms^^ 


n{R2r 
10^ cm"^ 


5'l.2mm 

MJy sr"^ 


N{H2T 
lO^^cm"^ 


-tkin 

K 


References 


L 183 


15:54:08.80 


-02:52:44.0 


(0,0) 


2.37 


20 


18.3 


7.2 


7.0 


(1) 


L 1544 


05:04:16.90 


+25:10:47.0 


(0,0) 


7.20 


14 


17.2 


6.7 


7.0 


(2) 


Oph D 


16:28:30.40 


-24:18:29.0 


(0,0) 


3.35 


8.5 


14.6 


5.7 


6.0 


(3) 


L 1517B 


04:55:18.80 


+30:38:03.8 


(-10,-20) 


5.87 


2.2 


7.1 


2.8 


10.0 


(4) 


L 310 


18:07:11.90 


-18:21:35.0 


(30, 80) 


6.70 


0.9 


7.5 


2.9 


9.0 


(5) 



Offsets of the eontinuum peak with respcet to the reference position. 
^ Systemic LSR velocity. 
^ Peak particle number density. 

Column density computed assuming Tdust — 8 K and k^, — O.Olcm^g"^. The values increase by 25% if Tdust — 7 K. 

*^ Gas kinetic temperature fr om t he literature. 

fP I Pagani et al. I (12003.120071). (2 ) ITafalla eTall l l2002l '). ICrapsi eFal] | [2007l ) . O HWard-Thompson et al.H1999l) . IHariu eral] 120081) (4) ITafalla et al.l 
j2004 l f5) IBacmann et al.l j2003) 



containing species in prestellar cores. Section [H] summarizes 
the model that has been adopted of the early stages of the 
collapse of the representative prestellar core L 1544. The 
fractional chemical abundances predicted by the model are 
presented and the corresponding column density profiles 
are compared with the observations. Finally, in Section [71 
we make our concluding remarks. 

2. Observational procedures and data reduction 

The observations were performed at the IRAM 30 m tele- 
scope in January 2008. The data have b een reduced and 
anlay zed using the CLASS90 software (jHilv-Blant et al.l 
20051) . The instrumenta l setup was identical to that used by 
Hilv-Blant et al.l (120081 ) : frequency-switching spectra, with 
a frequency-throw of 7.8 MHz, were recorded by the VESPA 
facility, with 20 kHz spectral resolution and 20 to 80 MHz 
bandwidth. The half-power beam-width is calculated as 
HPBW["] = 2460/i/[GHz], that is 28" at 87 GHz and 22" 
at 113 GHz. Cross-like patterns with 20" spacing were ob- 
served towards each source. The crosses were centered on 
the dust emission peak, as determined from published con- 
tinuum maps; all offsets quoted in the present paper refer 
to the central positions listed in Table [1] Table IC.ll sum- 
marizes the lines observed for each object. The amplitude 
calibration was checked every 10 min, the pointing every 
hour, and the focus every two hours, typically. Instrumental 
spectral effects were compensated by subtracting polyno- 
mial baselines from each spectrum before folding. More 
details on the data reduction procedures can be found in 
Appendix |B] All results (unless explicitly stated) have been 
translated from the antenna temperature scale (T^) to the 
main-beam temperature scale Tmb ~ x Fes/BcS, with 
the values of Bcs listed in Table IC.ll Fes is the forward 
efficiency, and Soft is the beam efficiency. 

3. Observational results 

All objects in our sample are pre-stellar, in the sense that 
none shows signposts of embedded stellar objects. Their 
peak particle number densities span more than an order of 
magnitude (see Table [J): the peak densities decrease from 
approximately 10^ cm^^ (L 183, L 1544) to 0.9 xlO^ cm^^ 
in L 310. All the objects have been extensively observed, 
both in their lines and continuum (see Fig. [T]) . Our obser- 
vations focused on the nitrogen-bearing species CN, HCN, 



"CN, Hi^CN, HNi^C, and HC^^N. AU these molecules 
present hyperfine structure (HFS). However, in the cases 
of HN^'^C and HC^'''N, the hyperfine structures were not 
resolved in our 20 kHz resolution spectra. All transitions 
are in the 3 mm band; the CN(2 — 1) line was observed in 
parallel at 1.3 mm (see Table ETTj) . 

3.1. Line properties 

All lines were observed successfully towards the three most 
centrally peaked cores, L 183, L 1544, and Oph D; the emis- 
sion lines are shown in Fig. [2] For CN and N2H+(1 — 0), 
only the weakest HFS components (at 113520.4315 MHz 
and 93176.2650 MHz) are shown. For each of the other 
lines, the strongest HFS component is considered instead: 
108780.2010 MHz, 86340.1840 MHz, 87090.8590 MHz, and 
86054.9664 MHz for "CN(1 - 0), H"CN, HN"C and 
HC15N(1 - 0) respectively. The HN^^C hyperfine structure 
is not fully resolved and the two hyperfine transitions of 
HC^^N are not resolved. The ^^CN line is detected in the 
four densest objects, L 183, L 1544, Oph D and L 1517B. 

Tables ID.lllD.5l give the properties of all lines towards 
all the positions in each source. The integrated intensi- 
ties, W, were derived from Gaussian fitting of a given HFS 
component (see above) . Several Gaussian components were 
fitted in some cases, e.g. for the known hyperfine struc- 
ture of HN^^C, and also in the obvious cases of multiple- 
component line profiles (L 1544). In such cases, W is the 
sum of the integrated intensities of each velocity compo- 
nent. The peak temperature Tmb is the maximum inten- 
sity over the line. Given that the lines are, in general, 
not Gaussian, the linewidth is estimated as the equivalent 
width. Aw = W/Tj^h and the statistical uncertainty is ob- 
tained by propagating the errors. For non-detections, upper 
limits on the integrated intensity were obtained by fitting a 
Gaussian at a fixed position. Upper limits on the intensity 
are at the 3a level while those on the integrated intensity 
are at the 5a level. 

Towards L 1544, all resolved lines show two clear peaks, 
with a dip centred at 7.20 kms~^. These t wo peaks are 
seen in several tracers including H^'^CO"'" bv iHirota et al.l 
(i20l)3h who concluded that there are two distinct veloc - 
ity components along the line of sight (jTafalla et al.|[l"998h . 
Owing to their double-peak profiles, lines towards L 1544 



Hily-Blant et at Nitrogen chemistry in starless cores 



3 




ax [seel ilx [seel 

Fig. 1. Continuum emission (MJy sr"-'^) at 1.2mm with the locations of the line integrations (crosses). The HPBW at 
1.2 mm and at the 3 mm frequencies discussed in this paper arc indicated. A linear scale of 0.05 pc is also shown, 
assuming a distance of 15 pc for all cores except L 183 ( 110 pc). Cont i nuuni maps for L 1544, Qp h D, L 183, L 1517B 
and L 310 are taken from IWard-Thompson et all (|l999t) . iPaeani et all ([2001), iTafalla et all (|2004 ) and iBacmann et all 
(|2000f) respectively. 



have the largest integrated intensities of all the lines that 
we observed. 

The comparison of the ^■^CN lincwidths shows that the 
lines towards L 183 are the narrowest with full widths at 
half maximum (FWHM) « 0.20 kms-\ The H"CN hue 
in this source exhibits a blue wing and the profile can be 
well fitted by two Gaussian components with FWHM of 
0.38 and 0.85 kms^^; this blue wing is not evident in any 
other tracer. Towards Oph D and L 1517B, the lincwidth 
is larger by a factor of 2 to 3, although the comparison 
with L 1544 is rendered difficult by the double-peak line 
profiles. In several tracers, e.g. H^'^CN, L 310 displays the 
largest linewidth (< 0.8 kms~^) but small integrated in- 
tensities. In all the sources, HC^^N has been detected, and 
the properties of the line, averaged over all offset positions, 
are listed in Table |21 The peak and integrated intensities 
decrease as the peak density decreases. Again, L 1544 seems 
to be an exception, owing to the double-peak line profile. 
The FWHM are comparable (w 0.4 kms~^) for all sources; 
but, once again, the FWHM is significantly larger (by a 
factor 2) in L 310 than in the other sources. 

3.2. Line ratios 

The ratios of total integrated intensities Wtot for some line 
combinations in each source, are shown in Fig. IC.ll Under 
the optically thin assumption, these ratios refiect the rel- 
ative abundances. The ratios CN/HCN or i^CN/H^CN 



are constant to within a factor of 2 across all the cores 
and vary between about 0.5 to 5 from source to source. 
Towards L 1517B, the ^^CN/Hi^CN ratio appears to de- 
crease towards the centre. Significant also is the fact that 
the H^'^CN/HN^'^C is constant and of similar magnitude 
(0.2-0.8) in all sources, independent of the central density. 

Most of the lines that have been observed are split by 
the hyperfine interaction, and the relative intensities of the 
hyperfine components can be used as a measure of opti- 
cal depth. It is generally assumed that the level popula- 
tions of hyperfine states are in LTE and hence proportional 
to the statistical weights, within a given rotational level. 
However, it has been known for some time that this as- 
s umption is often invali d (sec, for example, the discussion 
of lWalmslev et anil982l) , and this is confirmed by our data. 
When the populations are in LTE, one expects the satel- 
lite line intensity ratios to lie between the ratios of the line 
strengths, in the optically thin limit, and unity for high 
optical depths. As may be seen from Fig. [3l this is usu- 
ally but not always the case. For example, it is clear that 
the CN(1 — 0) ratios towards L 183 are inconsistent with 
this expectation, whereas the ratios observed towards other 
sources suggest high optical depths and are broadly con- 
sistent with equal excitation temperatures in the different 
components. The observed ^"^CN ratios show clear signs of 
deviations from LTE, but the effects are much less dras- 
tic than in the more abundant isotopomer, and we suspect 
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Table 2. Properties of the average HC^^N(1 — 0) profiles and average fractional abundance towards the observed cores 
(main beam temperature scale). 



Source 




(W) 




(FWHM) 


(iV(HC'-'N)) 


(iV(H2))20" " 


(iV)/(2iV(H2)) " 




mK 


mK kms^^ 


kms^^ 


kms^^ 


xlO^°cm-2 


xlO^^ cm-^ 


xlO"^^ 


L 183 


125 ± 15 


47 ± 4 


2.49 


0.35 ± 0.04 


8.9±0.8 


5.6 


0.8 


L 1544" 


170 ± 12 


85 ± 5 


7.20 


0.48 ± 0.05 


16.1±1.0 


4.3 


1.9 


Oph D 


102 ± 15 


47 ± 6 


3.44 


0.43 ± 0.08 


8.9±1.1 


5.6 


0.8 


L 1517B 


82 ± 15 


29 ± 3 


5.69 


0.33 ± 0.04 


5.5±0.6 


2.0 


1.4 


L 310 


20 ± 10 


21 ± 3 


6.69 


0.98 ± 0.16 


4.0±0.6 


0.8 


2.5 



" Using a two— compon ent Gaussian fit. Tiie FWHM is taken to be the equivalent linewidtli, in tliis case; vq — 7.2 kms ^ is adopted from 
IHilv-Blant et al.l Il2008fl . 

Column density of H2 derived from the dust emission smoothed to 20" , assuming T^ugt — 8 K (see TablefTl. The total {statistical and systematic) 
uncertainty is 30%. 
^ The uncertainty is typically 40%. 



that optical depths are low. In the case of H^^CN, depar- 
tures from LTE appear to be minor. 

The reasons for departures from LTE such as seen 
in Fig. 12] are presently unknown and need to be estab- 
lished. Such an investigation would require calculations 
of the collisional rate coe fficients, analogous to those of 
iMonteiro fc Stutzkil (jl986f ). as well as a treatment of the 
radiative transfer; this is beyond the scope of the current 
study. For the present, we use low abundance isotopomers 
such as ^'^CN and H^^CN to trace abundance gradients, 
neglecting collisional excitation and the possibility of frac- 
tionation of the ^^C and ^''N isotopomers. 



4. Column densities and abundance ratios 

We determine column densities using the standard formal- 
ism described in Appendix [C] (Eq. IC.ll and IC.2[) and as- 
sume the local (solar neighbourhood) ^^C :^^C= 68 ratio. 
It is instructive to consider also the abundance variations 
from source to source. Converting column densities into rel- 
ative abundances requires the molecular hydrogen column 
density, iYuj, which we have determined (indirectly) from 
measurements of the dust e mission, using bolometer maps 
available in t he lit erature (" Ward- Thomps on ct al l Il999t 
Pagani et ahl 120031: iTafalla eTal.l l2004t iBacmann et al.l 
2000ir and smoothing where necessary to a 20" beam. We 
assumed a dust temperature of 8 K and a 1.3 mm absorp- 
tion coefHcient of k = 0.01 cm^g-i (sec HWFP08). The 
results are shown in Fig. [5] for L1544, L 183, Oph D, and 
L 1517B. 

We remark first that, towards L 1544 , the column den- 
sities of ^^CN, H^^CN, and HN^^C are ah roughly pro- 
portional to the hydrogen column density, as inferred from 
dust emission (sec Fig. |4|); this has been noted already by 
HWPF08 for the case of "CN. The derived abundances 
do not change appreciably towards the dust peak, in spite 
of a variation of almost an order of magnitude in the hy- 
drogen column density. Thus, in this source, and with the 
current resolution, the CN-containing species do not ap- 
pear to be significantly depleted at high densities. On the 
other hand, the abundances tend to increase to the NW of 
the dust emission peak (see the SE-NW cut); we assume 
that this is related to asymmetry of the source. It is inter- 
esting that N2H"'" behaves in similar fashion. Important for 
the later discussion is the fact that the abundance ratios 



HN^^CiHi^CN and i^CNiH^^CN are approximately equal 
to 2 (with variations of up to a factor of 2); we assume 
that this reflects the ratios HNC:HCN and CN:HCN, re- 
spectively. 

However. L 1544 is not typical. Towards L 183, for exam- 
ple (sec the EW cuts in Fig.g]), the peak H^^CN and HN"C 
column densities are offset to the east, relative to the dust 
emission, whereas N2H+ and ^'^CN appear to follow the 
dust emission. The situation is similar in Oph D although 
we did not observe ^'^CN in this source. Towards L 1517B, 
HN^'^C correlates reasonably well with A^(H2) but this is 
not the case of ^"^CN nor, probably, of H^'^CN. Bearing in 
mind the inaccuracy of the abundance determinations, and 
the possibility of ^'^C fractionation, we conclude conserva- 
tively that there is no evidence for an order of magnitude 
variation in the CNtHCN nor the HNCiHCN abundance 
ratios between the dust emission peak and offset positions. 

We give also in Table [3] our estimates of the fractional 
abundances, relative to H, of CN, HCN, and HNC towards 
the dust emission peaks of our sample of sources; these 
abundances have been derived assuming the local value of 
68 for the ^^C:^^C ratio. All the relative abundances are 
of order 10~^, with CN and HNC being more abundant 
than HCN by a factor o f about 2. This va l ue is close to 
the ratio determined by llrvine fc Schloerbl (|l984f ) toward 
TMC-1. Values larger than 1 for the HNC:HCN abun- 
da nce ratio were found towards a sample of 19 dark clouds 
bv iHirota et al.l (jl998f ) with an average ratio of 2.1 ± 1.2. 
We do not see indications of significant abundance differ- 
ences between cores of high central density (L 1544 and 
L 183) and cores of lower central density (L 1517B, Oph D). 
Whilst the complexities of the source structure and of ra- 
diation transfer prevent our establishing the existence of 
small abundance differences, we may conclude that there 
remain appreciable amounts of CN, HCN, and HNC at den- 
sities above the typical density (3x10^ cm~'^) at which CO 
depletes on to grains. Not surprisingly, this effect is seen 
most readily in sources of high central column density, like 
L 1544, in which emission from the low density envelope is 
less important. 



5. Chemical considerations 

In this Section, we seek to update and extend previous 
studies of the interstellar chemistry of N-containing species 
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Fig. 2. Comparison of the line profiles of different tracers; the spectra have been shifted vertically for clarity. For CN(1— 0), 
(2 - 1), and N2H+(1 - 0), the weakest HFS components (at 113520.4315, 226887.3520 and 93176.2650 MHz respectively) 
are shown. For each of the other lines, the strongest HFS component is considered: 108780.2010, 86340.1840, 87090.8590, 
and 86054.9664 MHz for ^^Cn{l - 0), H^^CN, HN^C and HCi'^N(l - 0) respectively. The N2H+ spectrum is taken from 
A07. The spectra are for zero offset (cf. Table [B- Towards L 1517B, the CN(1 - 0) spectrum is replaced by CN(2 - 1) 
and the N2 H~'" spectrum is tak en from iTafalla et ahl (j2002t ). Towards Oph D, the N2H+ spectrum at offset (0", 0") is 
taken from lCrapsi et all ()2005() . 
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Table 3. Fractional abundances of CN, HCN, HNC and N2H+ relative to H towards the dust emission peaks of our 
source sample. Tables [DJIID.51 give the fractional abundances toward all the observed positions. 



Source 


[CN] 


[HCN] 


[HNC] 


[N2H+] 


[NO]' 


' CN:HCN 


HNC:HCN 




xlO-3 


xlO-« 


xlO"^ 


xlO-^° 


xlO" 


9 




L 183 


0.40 


0.17 


0.45 


0.78 


10.0 


2.4 


2.6 


L 1544 


1.11 


0.56 


1.12 


1.20 


4.0 


2.0 


2.0 


Oph D 


0.37 


0.23 


0.38 






1.6 


1.7 


L 1517B 


1.12 


0.48 


1.50 






2.3 


3.1 


L 310 


4.40 


1.50 


2.50 






2.9 


1.6 


An isotopic 


ratio ^^C : 


^^C = 68 was 


assume 


d when using 


the 


isotopologues . 





A global uncertainty (including statistical and systematic errors) of 40% is assumed. 
The abundance of NO is taken from A07. 



CN(l-O) 






1.1349 10 1.135 10 1.1351 10 1.1352 10 



Rest frequency [MHz] 



CN(l-O) 



1.0375 10° 1.0B735 10° 1.0579 10° 1.08795 10° 10' 



Rest frequency [MHz] 



H^^CN(l-O) 



B.6339 10 8.634 10 8.6341 10 3.6342 10 



l^est frequency [MHz] 



HCN(l-O) 



L310 

8.8631 10'' 8.8632 lO'' 8.8633 10* 8.8634 10* 



Rest frequency [MHz] 



Fig. 3. Relative integrated intensity of each HFS component for several species, at the central position for each source. 
In each panel, the dashed line indicates the relative intensities in LTE for optically thin emission. The thick line shows 
the optically thick limit. The abscissa is the rest line frequency. 



iPineau des Forets et all (|1990D : ISchilke et all (|1992D . with 
a view to providing a framework for the interpretation of 
our observations of prestellar cores. We shall show that it 
is possible to derive a simple expression for the CN:HCN 
abundance ratio, in particular, by identifying the principal 
reactions involved in the formation and destruction of these 
species. 

5.1. Main chemical reactions 

The fractional abundances of gas-phase species in prestellar 
cores are determined by: 

— the initial composition of the molecular gas which un- 
dergoes gravitational collapse; 

— variations of the density with time; 

— the rates of gas-phase reactions at the low temperatures 
(T « 10 K) of the cores; 

— the rates of adsorption of the constituents of the gas on 
to grains. 

The assumption of steady state, when computing the 
initial abundances, is not crucial when dealing with species 



which arc produced in ion-neutral reactions. The timescales 
associated with ion-neutral reactions are much smaller, in 
general, than dynamical timescales, notably the gravita- 
tional free-fall time. Consequently, the evolution of the 
fractional abundances of species which are formed in such 
reactions soon becomes independent of the initial values. 
However, it is believed that the timescales for producing N- 
containing species, such as NH3 and CN, are determined by 
slower neutral-neutral reactions. It follows that their frac- 
tional abundances in the subsequent gravitational collapse 
may depend significantly on the initial values. 

The conversion of atomic into molecular nitrogen in the 
gas phase is believed to occur in the reactions 



and 



The reactions 



N-fOH - 


NO + H 


(1) 


N-hNO - 


->N2 + 


(2) 


N + CH - 


CN + H 


(3) 


N + CN - 


-^N2 + C 


(4) 


C + NO - 


CN + 


(5) 
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Fig. 4. Derived column densities towards each source. From top to bottom: L 1544, L 183, Oph D, and L 1517B. The 
column densities are plotted as a function of the distance from the dust emission peak, along both cuts. Also plotted 
is the H2 column density (grey histogram, right scale), as derived from the dust emission (assuming T^ust = 8 K and 
K = 0.01 cm^ g~^). The right panels show the corresponding abundances. 



C + NO — ^CO + N (6) 
can also destroy NO, producing CN in the case of ([5]), and 

CN + O — > CO + N (7) 

can destroy CN. 

From the above, we see that NO forms from the reaction 
of N with OH, whereas CN forms from N and CH. It follows 
that the ratio of carbon to oxygen in the gas phase is a 
factor determining the relative abundance of NO and CN. 
The NO:CN abundance ratio will be lower in gas which is 
depleted of oxygen, either because of an intrinsically high 



C:0 elemental abundance ratio, or due to the diiferential 
freeze-out of O and C on to the grains, where the oxygen 
is incorporated mainly as water ice. 

Once N2 has formed, in ([5]) and N2H+ is produced 
in the protonation reaction 

N2 + H+ — ^ N2H+ + H2 (8) 

Dissociative ionization of N2 by He+ results in the produc- 
tion of N+: 

N2 + Hc+ — > N+ + N + He (9) 

Whilst the reaction of N+ with para-H2 (in its ground 
rotational state) is endothcrmic, by approximately 168 K, 
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Fig. 5. Same as Fig. S] for tlie derived fractional abundances. 



its reaction with ortlio-H2 is sliRhtly ex o therm ic and occurs 
even at low temperatures iLe BourlotI ()1991h . Subsequent 
hydrogenation reactions with H2 lead to NH^, which can 
dissociatively recombinc to produce NH3. Thus, N2 is a 
progenitor of both N2H+ and NH3, whilst NO and CN are 
intermediaries in the formation of N2 . 

HCN and HNC are produced principally in the reactions 



NHn 



■N 
C 



HCN 
HNC 



H 
H 



(10) 
(11) 



in which the products have so much excess energy that 
rapid isomerization is expected to yield practically equal 



amounts of HCN and HNC (jHerbst et al.ll2000t) . As CH and 
CH2 are produced through the dissociative recombination 
of hydrocarbon ions, notably CH;^, they are expected to 
coexist in the medium. It follows that CN, HCN and HNC 
should coexist also. HNC converts to HCN in the reaction 



HNC -f H^ 



HCN + H^ 



(12) 



The reverse reaction is endothermic and proceeds at a neg- 
ligible rate at low temperatures. HCN is destroyed princi- 
pally in the charge transfer reaction with H"*" and in the 
proton transfer reaction with Hj^ 



HCN -I- H^ 



HCN^ 



H 



(13) 
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Fig. 6. The principal reactions involved in the low- 
temperature chemistry of nitrogen-containing species, and 
specifically of CN, HCN and HNC. 

HCN + H+ — > H2CN+ + H2 (14) 

HCN+ reacts rapidly with II2, producing H2CN+, which 
dissociatively recombines with electrons, producing HCN 
and HNC; but there exists a branch to CN 

HjCN + e" — ^CN + H2 (15) 

for which the branching ratio /cn was taken equal to 1/3, 
and this represents a true destruction channel for HCN (and 
formation channel for CN). 

Because CN but not HCN (nor HNC) is destroyed by 
O, the abundance ratio CN:HCN increases as the C:0 ratio 
increases. The chemical network is summarized in Fig.l 

5.2. Simplified analysis of the chemistry 

We show now that an analysis of the principal reactions 
leading to the formation and destruction of CN and HCN 
leads to a simple analytical formula for the CN:HCN abun- 
dance ratio. When reactions (jSj and (|4|) determine the 
abundance of CN, 

n(CN) = n(CH);^/g]= 0.75n(CH) 

with the adopted values of the rate coefficients for these 
reactions at T = 10 K. 

We have seen already that HCN is formed in reaction 
pU)) . with a rate coefficient ^[Jq]= 2.2 x 10~^^ cm'^ s~^, and 
it is removed principally by H+ and H;^, forming H2CN+ 
(reactions [13] and [13]) . It follows that 

A[i2p(CH2)n(N) = fcL?i(HCN)[n(H+) +n(H+)]/cN, 

where is the Langcvin rate coefficient which character- 
izes the reactions of HCN with H+ and Hj^. 



Under the above assumptions, the CN:HCN ratio is pro- 
portional to the CH:CH2 ratio: 

MCN^^^niCHl fcLn(e-) ^ n(CH) n(e-) 
n(HCN) ^7i(CH2)-^™^(N) ^ 7i(CH2) ?i(N) 

where wc have assumed that n{e~) = n(li~^) + n{}i^). 
CH has its chemical origin in the reaction of He+ with 

CO 

He+ + CO — >C+ + + He (16) 

for which the adopted value of the rate coefficient is = 
1.5 X 10^^ cm'^s^^ (essentially Langevin). As CO depletes 
on to grains, there is competition with reaction p6l) from 
the reactions of He^ with H2, 

Hc+ + H2 — > H+ + H + He (17) 

He+ + H2 — > H+ + He 

4.0x10-" cm^s-i (jSchauer et al.lll989t) . Thus, the 
rate of production of carbon ions is CHe'T-(He)fc, where 
Chc ~ 0.5Ch9 is the rate of cosmic ray ionization of He 
and /c = ^(CO)/[/27Fi(H2) + ^^(CO)]. 

Most ot the C+ ions produced in reaction (|T6)) combine 
radiatively with H2 to form CH^ , which reacts rapidly with 
H2 , forming CHj^ which then recombines dissociatively with 
electrons, yielding CH 

CH+ + c^ — > CH + H2 (18) 

CH++C- — ^CH + H + H (19) 

but also C and CH2. The ratio feu = 0.3 is the fraction 
of the dissociative recombinations of CHjj" which form CH 
and is significant for the CN:HCN abundance ratio. CH and 
CH2 are destroyed by atomic O: 

CH + O — > COH (20) 

CH2 + O — >CO + H + H (21) 

and 

CH2 + O — >C0-|-H2 (22) 

forming CO, and hence the elemental C:0 ratio is, 
once again, a pertinent parameter for the CN:HCN 
ratio. So too are the rate coefficients for reac- 
tions (dOl), and ([21]). We adopted the (temperature- 
independent) values of these rate coefficients in the 
NIST chemical kinetics databas^B- All three reac- 
tions are rapid, with rate coefficients of the order of 
lO"^*' cm^ s^^. The number density of CH is predicted 
to be n(CH) = CHen(He)/c/cH/[^(N) + /^(O)], 
Similarly that of CH2 is predicted to be n(CH2) = 
CHcn-(He) fcfcR2/[I^Up'(^) + (to]+ te])"-(0)]- Assuming 
that tex] ^ ^15^ "^a^ deduce from the above equations 
that the CH:CH2 ratio is given by 

njCR) ^ /cH 2/^(0) + /^(N) 
«(CH2) /cH. ^(0) + /gp(N) 

It follows that n(CH)/n(CH2) depends on the ratio 
of the atomic N and O abundances and varies between 
extreme values of approximately 2, when n{0) ^ 't-(N), 
and 0.2, when 7i(N) ^ 7i(0). Model calculations (see 

^ http : / /kinetics . nist . gov/kinetics/ 
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Section r6.2p confirmed that the above approximation repro- 
duces the calculated CN:HCN abundance ratio to within a 
factor of approximately 2 over most of the density range 
cm ^. In the following Section we con- 
sider the issue of the atomic nitrogen abundance in prestel- 
lar cores in more detail. 



5.3. The abundance of atomic nitrogen in cores 

Atomic abundances in prestellar cores are notoriously diffi- 
cult to determine. Although the atomic fine structure tran- 
sitions are, in principle, observable, it is difficult, in prac- 
tice, to distinguish a component corresponding to dense, 
cold molecular material from emission arising from low den- 
sity, hotter layers along the line of sight. The emission from 
photon dominated regions (PDRs), for example, tends to 
be stronger than that from cores. 

The analysis in Section suggests that the CH:CH2 
ratio is dependent on the abundances of atomic oxygen and 
nitrogen. Combining the expressions for the CH:CH2 and 
the CN:HCN ratios, there follows the inequality 



x(N) = ?i(N)/n(H2) < 100 



n(CH) x{e) 
7i(CH2) i?(CN : HCN)^ 



where x(e) denotes ri(e)/n(H2) and i?(CN : HCN) is the 
measured value of the CN:HCN ratio. The inequality arises 
from neglecting the destruction of CN by oxygen in reaction 
([7]). We a dopt n(CH)/n(CH2) ^ 1 and estimate x{e) from 
Fig. 5 of IWalmslev et al.l (I2004D as 5 x lO-^n^o ^ where 
ns is n(H2) in units of 10^ cm~'^. From Table [31 we infer a 
typical value of 2 for i?(CN : HCN) and conclude that the 
fractional atomic N abundance 

"(N)/nH < 1.25x10"^ 

This value is considerably less than the cosmic nitro- 
gen a bundance (n(N)/nH = 7.9 xlO~^; I Anders fc Grevessd 
llQSgT) and implies that there is only a small fraction of el- 
emental nitrogen in gaseous atomic form. Although some 
approximations were made when deriving this value (such 
as n(CH)/n(CH2) « 1), our more precise and complete 
chemical modelling has confirmed that this estimate is es- 
sentially correct as long as the total gaseous nitrogen abun- 
dance is n (N)/nH < 10"^ (see Fig. [3). It follows that most 
of the elemental nitrogen must be in the form of gaseous 
N2 or N-bearing solid compounds {e.g. NH3 or N2 ices). 
Indeed, for n(N)/nH < 10~^, steady-state models suggest 
that most of the nitrogen is in ices (see Section |6] as well as 
the discussion of Maret et al. 2006). 

An analogous inequality, in terms of the measured 
CN:HCN ratio, can be derived for the fractional abundance 
of atomic oxygen, a:;(0). However, it is a much weaker con- 
straint than the limit on x(N). Other observables, such as 
NO, provide stronger constraints on the fractional abun- 
dance of atomic oxygen (see A07) . 



6. Models 

6.1. Steady state 

The timescale for the nitrogen chemistry to reach steady 
state is known to be large, relative to the free-fall time in a 
prestellar cor e, owing to the slow conversion of N to N2 (see, 
for example. iFlower et al.l [20061 ). Our model calculations 
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Fig. 7. Steady-state fractional abundances of nitrogen- 
containing species for a density njj = 10* cm""^, a ki- 
netic temperature T = 10 K, and a cosmic ray ioniza- 
tion rate of C = 10~^^ s~^. Upper panel: C:O=0.97 is held 
constant and the gas-phase nitrogen fractional abundance 
[N]gas = n(N)/nH is varied from 10~^ to 6.4x10"^. Lower 
panel: [N]gas = 10~^ is held constant and the gas-phase 
C:0 ratio is varied by changing [Ojgas- 



show that, with a cosmic ray ionization rate C = 10^^^ s~^, 
this timescale is of the order of 10^ yr. As a consequence, 
the results of the time-dependent models of gravitational 
collapse depend on the initial composition which is adopted 
and the rate of dynamical evolution. 

Following the discussion in Section [5l it is nevertheless 
instructive to examine the results of steady-state calcula- 
tions (that is to say time independent and only in the gas 
phase), as functions of the fractions of elemental oxygen 
and nitrogen in the gas phase. In this way, an impression 
may be obtained of the dependence of the observables on 
the degrees of depletion, without the complications of the 
time dependence, which is considered in Section [621 

In the upper panel of Fig. [71 we present the fractional 
abundances of nitrogen-containing species in steady state 
for a density tt-h — 10^ cm~'^, a kinetic temperature T = 



10 K, and a cosmic ray ionization rate of C = 10' 
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For reaction (|7]), we adopted a temperature-independent 
rate coefficienlQ (4 x 10"^"'^ cm^ s""'^) although we note that 
there is some theoretical evidence th at the rate of this reac- 
tion m ay decrease with temperature (jAndersson fc Markovil 
|2003[ ). The fraction of elemental nitrogen in the gas phase 
varies from 0.017 to 1. We assume implicitly that the 'miss- 
ing' nitrogen is in solid form. Following A07, we adopt a 
relative abundance of elemental carbon to oxygen in the 
gas phase C : O = 0.97. In the lower panel, the fractional 
abundance of elemental nitrogen in the gas phase is held 
constant at rt(N)/nH = 10~^ and the C:0 ratio is varied 
by changing the fractional abundance of oxygen in the gas 
phase [0]gas. 

We see from Fig. [7] that, in steady state, there tends to 
be somewhat more molecular than atomic nitrogen in the 
gas phase. Species such as NH3 and N2H+ have abundances 
which are roughly proportional to N2 . The HCN abundance 
is relatively insensitive to changes in the gas phase C:0 
ratio but follows the gas-phase nitrogen abundance. On the 
other hand, CN and NO are sensitive to the C:0 ratio. The 
net effect is that CN:HCN increases with C:0 and decreases 
with the fraction of nitrogen in the gas phase. 

6.2. Gravitational collapse 

A more satisfactory approach to comparing observational 
results with models is through a simulation of a gravita- 
tional collapse. In Fig. |S1 wc compare the observed values 
of the HCN abundance and the CNiHCN ratio with the pre- 
dictions of models in which the density and the chemistry 
evolve following free-fall gravitational collapse. All neutral 
species are assumed to adsorb on to dust grains (of ra- 
dius 0.5 ^m) with a sticking coe fficient of unity an d are 
desorbed by cosmic ray impacts (iFlower et al.ll2006l ). We 
assume that, initially, the chemical composition of the gas 
has attained steady state at a density nn = 10* cm~^. We 
make various assumptions concerning the amount of nitro- 
gen initially in the gas phase (or, equivalently, the fraction 
which is initially in the form of nitrogen-containing ices on 
grain surfaces). 

The fraction of elemental nitrogen in the ambient molec- 
ular medium which is in solid form is poorly known. There 
is evidence for ammonia ice in spectral profiles observed 
towards some youn g YSOs, with per haps 15% of the abun- 
dance of water ice (jGibb et al.ll2000t ). but no such evidence 
exists towards background stars; there is perhaps a substan- 
tial fraction of the nitrogen in the form of N2 ice also. In 
prestellar cores, the abundance of N2II+ places lower limits 
on the amount of gas-phase nitrogen, which we estimate 
conservatively to be about 10^®. Accordingly, we have var- 
ied the initial gas phase nitrogen abundance in our models 
in the range 10"^ < «(N)/nH < 6.4 x 10~^, where the upper 
limit co rresponds to the valu e observed in diffuse interstel- 
lar gas (jSofia fc Meverl[200lt) . Fig. [8] displays the fractional 
abundances computed in the course of the collapse, at den- 
sities TT-H = 10^ cm~'^ and nn = 10^ cm~^. We show also, 
for comparison, the initial (steady state) values. 

It may be seen from Fig.|5]that, at a density of 10^ cm~^ 
which is the more relevant value for the purpose of com- 
paring with observations, reasonable agreement is obtained 
only for initial gas-phase nitrogen abundances close to 10^® 

^ From the osu_03_2008 rates of Eric Herbst's group 
(http : //www. physics . ohio-state . edu/'eric) 



- in other words, close to the lower limit. Even so, the com- 
puted abundances do not fit well the observations of L 183 
and Oph D; but wc note that the density at the dust peak 
in L 183 approaches 10^ cm~^. Our model results are de- 
pendent also on the fraction of oxygen locked in ices, or, 
equivalently, on the initial gas-phase C:0 ratio. It is possi- 
ble that this ratio varies considerably from source to source, 
resulting in discrepancies when we compare observations 
with model predictions. 

Figure[5] (bottom right panel) shows the HCN:CO abun- 
dance ratio as a function of the CO fractional abundance. 
For a given initial abundance of gaseous nitrogen [N]gas 
(= n(N)/nH), the abundance ratio is followed along the 
collapse and values are shown at densities tt-h = 10^, 4x10^ 
and 10*^ cm-3. The differential freeze-out of HCN and CO 
is evident. In all these models, CO depletes by two orders 
of magnitude. The behaviour of HCN regarding depletion 
is different in that it depends on the initial [Njgas- For a 
large initial [N]gas = 6.4x10^'"', HCN depletes only a factor 
of 3 less than CO. However, at the other extreme value 
([N]gas = 10"^), HCN depletes 10 ti mes less thaii CO 
Observational values towards L 1544 (jCaselli et al.lll999t) 
and L 1517B (jTafalla et al.ll2002l) favour differential freeze- 
out between HCN and CO and thus low initial [N]gas. 

We conclude from Fig. |S] that the models fail to explain 
the observations. One possible reason for this failure is our 
neglect of line-of-sight effects in the models used to con- 
struct Fig. [51 The observed quantities are column densities, 
which arc integrals along the line of sight over a range of 
densities; our analysis neglects this effect. However, trial 
calculations for one source (L 1544; see Appendix [D|) sug- 
gest that including line-of-sight integration can reduce but 
not eliminate the discrepancies between model predictions 
and observations. Another possibility might be that the 
duration of the collapse is longer that the free-fall time. 
However, if this time is significantly increased, the abun- 
danc e of gaseous CO dro ps too rapidly with increasing den- 
sity (jFlower et al.ll2005h . More important may be errors in 
the rate coefficients that we have used for some of the key 
reactions, discussed in Section [SJ It is clear, for example, 
that our predictions relating to CN are sensitive to the rates 
of reactions ([3]), (|3]), and pUj) at temperatures of the order 
of 10 K. Further progress in this field will require reliable 
determinations of the rate coefficients of these reactions at 
low temperatures. 

We find that the values of the CN:HCN ratio observed in 
prestellar cores indicate that the fraction of nitrogen in the 
gas phase is likely to be considerably lower than the diffuse- 
gas value of6.4xl0~^. Nitrogen (like oxygen) may deplete 
on to grain surfaces at relatively low densities. Confirmation 
will require the identification of nitrogen-containing ices 
and estimates of their relative ab undances. A r a ther s imilar 
conclusion has been reached by iMaret et al.l (j2006| ) in a 
study of B68. 

We finally note that our observations show HNCTlCNw 
2, whe reas the exhaustive theoretical study of lHerbst et al.l 
()200GD predicted HNC:HCN « 1. It is possible that en- 
hanced line trapping in HNC, relative to HCN, results in 
our deriving an anomalously high HNC:HCN abundance 
ratio; but it is unlikely that this effect can explain fully the 
discrepancies with the model predictions. Maybe more rel- 
evant is the apparent correlation between the HNC:HCN 
ratio and freeze-o ut, as suggested by the results from 
iHirota et al.l ([l99l) who show that the largest (resp. small- 
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Fig. 8. The fractional abundances computed by the gravitational collapse model. We show in panel (a) results for the 
HCN abundance as a function of the CN:HCN ratio. In panel (b) we show the HCN abundance as a function of the N2H^ 
abundance. In panel (c) the HCN abundance as a function of the NO abundance. The initial C:0 ratio is assumed to be 
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L 1517B are taken from lCrapsi et all (|2005[ ). In panel (d) we show the HCN:CO abundance ratio as a function of the CO 
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est) ratio is observed toward a strongly depleted core (resp. 
undeplctcd) . 



7. Concluding remarks 

We have studied the behaviour of nitrogen-containing 
species, principally CN, HCN, and HNC, in the pre- 



protostcUar cores L 183, L 1544, Oph D, L 1517B, and 
L 310. Our main conclusions are as follows. 

— We observe that CN, HCN, and HNC remain present 
in the gas phase at densities above the typical density 
(3x10"^ cm~^) at which CO depletes on to grains. 

- The CN:HCN and HNC:HCN ratios are larger than 
unity in all objects and do not vary much within in 
each core. Whilst the differential freeze-out of CN and 
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CO can be understood, the approximate constancy of 
the CN:HCN ratio cannot. 
— The CN:HCN ratio puts upper hniits on the abundance 
of atomic nitrogen in the gas phase, and the NO:HCN 
ratio constrains the C:0 ratio. Though uncertain, the 
comparison between observations and models indicates 
that most of the nitrogen is locked into ices, even at 
densities probably as low as 10"* cm~^. 

Our current knowledge of the chemistry of nitrogen- 
containing species is incompatible with the observed ra- 
tios. We recall that there exist large uncertainties in the 
rate coefficients, at low temperatures, for the key neutral- 
neutral reactions, discussed in Section [SJ including those 
with atomic nitrogen. Essential to further progress in un- 
derstanding the chemistry of nitrogen-containing species 
in pre-protostellar cores are measurements, at low temper- 
atures, of at least some of these key reactions. 
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Appendix A: Spectroscopic data 



Table A.l. Hyperfine struc t ure in CN(2 ~ 1) {NJF 
N'J'F'); from lSkatrud et all (Il983h . 



N 


, J 


F - 


->iV' 


, J' , F' 




R.I." 


2 


,3/2 


1/2- 


1 


,1/2,1/2 


226663.685 


0.0494 


2 


,3/2 


3/2- 


1 


,1/2,1/2 


226679.341 


0.0617 


2 


,3/2 


1/2- 


1 


,1/2,3/2 


226616.520 


0.0062 


2 


,3/2 


3/2- 


1 


,1/2,3/2 


226632.176 


0.0494 


2 


,3/2 


5/2- 


1 


,1/2,3/2 


226659.543 


0.1667 


2 


,3/2 


1/2- 


1 


,3/2,1/2 


226287.393 


0.0062 


2 


,3/2 


3/2- 


1 


,3/2,1/2 


226303.049 


0.0049 


2 


,5/2 


3/2- 


1 


,3/2,1/2 


226875.896 


0.1000 


2 


,3/2 


1/2- 


-> 1 


,3/2,3/2 


226298.896 


0.0049 


2 


,3/2 


3/2- 


1 


,3/2,3/2 


226314.552 


0.0120 


2 


,3/2 


5/2- 


1 


,3/2,3/2 


226341.919 


0.0053 


2 


,5/2 


3/2- 


1 


,3/2,3/2 


226887.399 


0.0320 


2 


,5/2 


5/2- 


1 


,3/2,3/2 


226874.183 


0.1680 


2 


,3/2 


3/2- 


1 


,3/2,5/2 


226332.519 


0.0053 


2 


,3/2 


5/2- 


1 


,3/2,5/2 


226359.887 


0.0280 


2 


,5/2 


3/2- 


-> 1 


,3/2,5/2 


226905.366 


0.0013 


2 


,5/2 


5/2- 


1 


,3/2,5/2 


226892.151 


0.0320 


2 


,5/2 


7/2- 


1 


,3/2,5/2 


226874.764 


0.2667 



a; Rest frequency in MHz. 

b: Relative intensities, with their sum normalized to unity. 



Appendix B: Data reduction 

Data reduction was done with the CLASS90 software from 
the GILDAS program suitfl We summarize here the re- 
duction of the frequency-switched spectra obtained with 
the VESPA autocorrelator facility at the IRAM 30 m radio 
telescope. 

All spectra were corrected first from the instrumental 
spectral transfer function. The resulting spectra were then 
folded and averaged (each folded spectrum being weighted 
by the effective rms of the residuals after baseline subtrac- 
tion). A zero-order polynomial was fitted to the resulting 
spectrum to compute the final rms. Whenever platforming 
was present in the data, the spectrum was split into as many 
parts as needed, and each part was treated individually 
with a first-order polynomial in order to adjust the contin- 
uum level. The concatenated sub-parts were then treated 
as a single spectrum. This method proved to be robust. In 
some cases, the amplitude of ripples was large enough to re- 
quire special treatment. Ripples were subtracted by fitting 

^ Available at http://www.irain.fr/GILDAS 
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a sine wave to the spectrum, using an improved version of 
sine fitting, as compared with the defauh CLASS90 proce- 
dure. In the case of a spectrum presenting ripples, channels 
in the spectral windows (where there is presumably some 
line emission) were replaced by a sine wave, determined by 
the first-guess parameters (amplitude, period and phase). 
As the critical parameter in the sine wave minimization 
proved to be the period, the minimization was repeated 
for several values of the period. In all cases, this algorithm 
converged to an acceptable solution, as indicated by the 
residuals and inspection by eye. 



Appendix C: Column density derivation 

All column densities are derived assuming optically thin 
emission with levels populated in LTE at the excitation 
temperature Tex. Under these assumptions, the column 
density is directly proportional to the integrated flux in the 
line Wtot- From a transition Ju J\ between energy levels 
and E\ (corresponding to an energy Tq = {E^ — E\) / ks.) , 
one can compute the total column density of the molecule 
as 

3eo/i 1 ge^>/'=Bi^- Wtot 



A^t, 



27r2 Ju/i2 AJ^[1 - exp(- 



(C.l) 



in SI units, where AJ,y = J^{Tcx) — J}^{Thg), with J,y{T) = 
To/[exp(To/T) — 1], and fio is the dipole moment; eg is the 
permittivity of free space. Numerically, we obtain 



8x10 



12 



JuA^D AJ,[l-exp(-^)] 

(C.2) 

with Wtot in Kkms ^and A^tot in cm is the dipole 

moment in debye. We have introduced the conversion fac- 
tor, A'o, which depends on the molecular properties and 
the excitation temperature. It is to be noted that, in the 
case of resolved hyperfine structure, Wtot is the total inte- 
grated intensity of the hyperfine multiplet, i.e. of the rota- 
tional transition. Values of Nq arc given in Table IC.ll for 
the observed molecules and at an excitation temperature 
= 8 K (see HWFP08). 



Appendix D: Modelling of the prestellar core 
L 1544 

In order to interpret our observations of L 1544, we make 
use of the model of one-dimensional, free-fall gravitational 
collapse used in our previous study (A07). This model in- 
corporates dust grain coagulation and a time-dependent 
chemistry, including the reactions listed in Section[5]above, 
which are directly relevant to the present work. We assume 
a constant kinetic temperature, T = 10 K, and a cosmic 
ray ionization rate, ^ = 10 ~^'' s~^. Furth er information on 
the model may be found in lFlowe~ et al.l dlOOS). 

An important aspect of the interpretation is to connect, 
as realistically as possible, the abundance profiles (i.e. the 
number density of species X, n(X) as a function of the total 
density, tih = 2n(Il2) -I- n(H)), which are the output of 
the model, to the observed variations in L 1544 of column 
densities, A'^(X), with impact parameter r. In order to make 
this connection, we proceed as follows: 



we relate the gas density at r to the central density by 
means of the relation 

nH(r) =nH(0)/[l + (r/r2)"] 



iTafalla et all (|2002f ). where r is the offset from centre, 
r = 0, and r2 is the radial distance ov er which the 
densit y decreases to nH(0)/2. Following iTafalla et al.l 
(|2002[ ). we adopt = 20" (equivalent to 0.014 pc at 
the distance of L 1544) and a ~ 2.5; the central den- 
sity n(0) = 10^ cm~^ (A07), which is somewhat smaller 
than the value reported in Table [T] for this object but 
within the probable uncertainties of its determination. 

— Using the computed values of n(X) vs nn, we calculate 
the corresponding column density, A^(X), by integrat- 
ing along the line of sight for any given value of r in the 
adopted range < r < 120 arcsec. over which the den- 
sity nH(f) decreases from 10^ cm~'^ to 1.1x10^ cm"-^. 

— Finally, the column densities are convolved with a 
Gaussian profile with a (1/e) radius of 15", correspond- 
ing to a HPBW of 25", in order to simulate approxi- 
mately the IRAM 30 m telescope beam at 100 GHz. 

As will be seen below, the consequence of this procedure 
is a significant damping of the variations in the computed 
abundance profiles, owing partly to the effect of integrating 
along the line of sight and hence over a range of densities, 
and partly to the Gaussian beam averaging. 

We consider first the predictions of the chemical model, 
and specifically the abundances of nitrogen-containing 
species. We turn our attention then to the Gaussian-beam 
averaged column densities, and their comparison with the 
observations. 



D.l. Abundance profiles 



In Fig. ID. II are plotted the fractional abundance profiles of 
CN, HCN, NO and N2II+; note that the x-axis has been 
reversed in order to facilitate the comparison with later 
Figures, in which the x-coordinate is the offset from the 
centre, where the density of the medium is highest. We see 
from Fig. ID. II that, at low densities, the fractional abun- 
dance of HCN exceeds that of CN, by a factor which ap- 



proaches two orders of magnitude when nn 



10^ 



The fractional abundance of HCN decreases towards the 
maximum density of 10^ cm~'^, where n(HCN) « n(CN). 
This behaviour can be understood by reference to the dis- 
cussion in Section [5] above: CN is formed and destroyed in 
reactions ^ |4]) which involve atomic nitrogen; HCN, on 
the other hand, is formed in reaction pH)) with N but de- 
stroyed in reactions with H+ and Hj^ that ultimately lead to 
CN. Consequently, as the density of the medium increases, 
and neutral species begin to freeze on to the grains, the 
fractional abundance of HCN falls, whereas the fractional 
abundance of CN remains roughly constant until, finally, 
CN too freezes on to the grains. 

We have already seen in Section [2] that the ratio of the 
column densities of HCN and CN, which is of the order of 
1 in L 1544, shows little variation across this and the other 
sources in our sample. Thus, it seems unlikely that the frac- 
tional abundance profiles of HCN and CN seen in the upper 
panel of Fig. ID. II are compatible with the observed col- 
umn densities, given that the ratio n(HCN)/7i(CN) varies 
by approximately 2 orders of magnitude over the range of 
density 10'* < ?t-h < 10^ cm"'^; this tentative conclusion is 
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Fig. C.l. Ratio of the total integrated intensities Wtot of the molecules indicated in each panel, as a function of the 
offset from the dust emission peak. For each molecule, Wtot is derived from the values in the corresponding Table divided 
by the relative intensity of the HFS component, when resolved. Filled and open circles distinguish between the two cut 
directions. Open triangles are 5(j upper limits. From top to bottom, L 183, L 1544, Oph D, L 1517B, and L 310. 
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Table C.l. Conversion factor, TVo (in 10^^cni^^/(Kkms~^)), at an excitation temperature Tex = 8 K (see Eg. IC.2I 



Molecule 


Transition 


a 

U 


ov 


DC 


B 


A*o 


T> T d 
H.l. 


Q 


iVo/R.I. 








MHz 


1 —1 




MHz 


debye 










(1 


-0) 




u.uoz 


n 7/1 


ODDyo.4 ( u 






0.0 


zyo 


CN 


(2 


-1) 


226882.000 


0.053 


0.50 












HCN 


(1 


-0) 


88633.936 


0.066 


0.77 


44315.976 


2.9852 


0.1111 


4.1 


16.4 


"CN 


(1 


-0) 


108780.201 


0.063 


0.77 


54353.130 


1.45 


0.194 


3.4 


29.5 


H"CN 


(1 


-0) 


86342.251 


0.054 


0.75 


43170.127 


2.9852 


0.5556 


4.2 


3.4 


HN"C 


(1 


~0) 


87090.850 


0.068 


0.78 


45331.980 


3.05 


1.000 


4.0 


1.7 


HC^^N 


(1 


-0) 


86054.9664 


0.067 


0.78 


43027.648 


2.9852 


1.000 


4.2 


1.9 


N2H+ 


(1 


-0) 


93171.621 


0.068 


0.78 


46586.880 


3.40 


0.037 


3.9 


35.2 



^ Frequency of the HFS eomponent considered for which wc compute the fiux W. 
^ Spectral resolution in kms^^. 

^ Beam efficiency. The forward efficiency at 3 mm is i^^ff — 0.95, and 0.91 at 1.3 mm. 
Relative intensity of the HFS component. The total fiux in Ea. [CT2l is Wtat = W/R.I. 

20 _ hBJ(J+l) 

^ Partition function computed as Q — y ^ (2 J -i- 1) exp f^u^^x ^ where B is the rotational constant. 

J=o 



confirmed by the analysis in the following Section ID. 21 It 
appears that the observations are indicating that the na- 
tures of both the formation and the destruction processes 
are similar for both these species. This realization led us 
to consider the possibility that reaction (U) , which destroys 
CN, may have a small barrier, significant at the low tem- 
peratures relevant here (T < 10 K) but difficult to detect 
by measurements at higher temperatures. A similar situa- 
tion may obtain for the analogous reaction ^ of NO with 
N. Accordingly, we show, in the lower panel of Fig. ID. II 
the results which are obtained on introducing a barrier of 
25 K to both reaction ([5]) and reaction we note that 
the adopted barrier size is arbitrary, the only requirement 
being that it should be small. 

Small barriers can arise when the potential energy 
curves involved in the atom-molecule reaction exhibit 
(much larger) barriers for certain angles of approach but 
no barrier for others. In order to determine the thermal 
rate coefficient, the probability of the reaction must be av- 
eraged over the relative collision angle. If the rate coefficient 
is then fitted to an Arrhenius form, 

A:(r) =7(r/300)"exp(-/3/T), 

small, positive values of /S may be interpreted as 
the angle-averaged value of the reaction barrier; see 
lAndersson fc Markovil ()2003[ ). 

It is clear from Fig. ID. II that the small reaction barrier 
has the effect of enhancing the fractional abundance of CN 
and reducing the amplitude of the variation in the ratio 
n(HCN)/n(CN). We shall see in the following Section |D2] 
that this variation is damped further when the ratio of the 
corresponding Gaussian-beam averaged column densities is 
considered. 

The results in Fig. ID. II have been obtained assuming 
that the grain-sticking probability was unity for all species, 
and that the elemental abundance ratio C : O = 0.97, i.e. a 
marginally oxygen-rich medium. This value of the C:0 ra- 
tio was an outcome of the modelling by A07 of observations 
of NO and NjE"*" in L 1544. These authors investigated also 
the consequences of varying the values of the sticking coef- 
ficient for atomic C, N and O, and the initial value of the 
N:N2 abundance ratio. We recall that increasing the ele- 
mental C:0 abundance ratio further has the consequence of 



10 ' , HCN 




iijj (cm"') 

Fig. D.l. The number densities of CN, HCN, NO and 
N2II+, relative to nn = 2n(Il2) -I- n{}i), as predicted by 
a model of gravitational free-fall collapse. The lower panel 
illustrates the effects of introducing a barrier of 25 K to 
both reaction ^ and reaction Qj. 



reducing the HCNiCN ratio, thereby improving the agree- 
ment with the observations. On the other hand, the CN:NO 
ratio also rises, and the values of this ratio in Fig. ID. II 
where n(CN) > n(NO), already exceed the values of the 
corresponding column density ratio, observed in L 1544, 
where iV(CN) < A^(NO) by typically an order of magni- 
tude. It is possible that selective variations in the values of 
the sticking probability, or in the initial N:N2 abundance ra- 
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beam averaging is a flattening of the column density pro- 
files (Fig. ID.2|) . compared with the fractional abundance 
Comparing the two panels of Fig. ID. 21 



profiles (Fig. ID.1[) 
we see that the introduction of the small barriers to the 
reactions of CN and NO with N reduces substantially the 
Ar(HCN)/7V(CN) column density ratio. Although they do 
not attain the observed value, of the order of 1, the com- 
put ed va lues of 7V(HCN)/iV(CN) in the lower panel of 
Fig. ID.2l are clearly more compatible with the observations 
of L 1544 than are those in the upper panel. 



Fig. D.2. The column densities of CN, HCN, NO and 
N2H+, relative to H2, as predicted by the model described 
in Section [Dl The lower panel illustrates the effects of in- 
troducing a barrier of 25 K to both reaction 1^ and reac- 
tion dH). 



ticQ, might alleviate some of these discrepancies. However, 
whilst there remain such large uncertainties in the values of 
the rate coefficients for the key neutral-neutral reactions, 
discussed in Section[Sl it would perhaps be premature to in- 
vestigate further the consequences of modifying the values 
of other (and equally uncertain) parameters, in an attempt 
to improve the agreement between the models and the ob- 
servations. Our aim here is to point to the discrepancies and 
highlight the uncertainties; and it seems unlikely that fur- 
ther progress can be made until the rates of at least some of 
the key reactions have been measured at low temperatures. 



D.2. Column densities 

In Fig. ID. 21 are shown the computed column densities of 
CN, HCN, NO and N2H+, relative to the column density 
of H2. The fractional abundances of these species, relative 
to riH, derive from the models discussed in the previous 
Section ID. II 

At zero offset, the line of sight passes through re- 
gions with densities covering the entire range of the model, 
l.OxlO^cm-3 > nn > 1.1x10'* cm"^ Consequently, the 
column density ratio, iV(HCN)/iV(H2) > 7i(HCN)/71h, 
evaluated at the peak density nn = 10® cm~'^; further 
smoothing is introduced by the Gaussian-beam averag- 
ing. The overall effect of the line-of-sight and Gaussian- 

* The equilibrium value of the N:N2 ratio, adopted in the 
present models, is n(N)/n(N2) = 1.0, as compared with the 
(non-equilibrium) value of 18 adopted by A07. 
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Table D.l. Line properties (main beam temperature scale) and total colmnn densities derived towards L 183. Numbers 
in parentheses are powers of 10. 
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mb 
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Af(H2) 


JVt„t/2Af(H2)s 








mK 
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0.26 


53 
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— 10) 
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— 10) 
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3.2 
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0.22 


38 
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60 
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0.21 
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4.2 
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20 
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11 
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0.49 


0.34 
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17 
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Error bars are Ict, and upper limits on W and A^tot are at the ba level. A final calibration uncertainty on the column density of 10% has been 
adopted, unless smaller than Icr. We adopt a systematic uncertainty on the dust column density of 30% which reflects the uncertainties on T^ust 
and K,jy . 

" Strongest component at 108780.2010 MHz, with R.I. ^ 0.194. 

^ Strongest component at 86340.1840 MHz, with R.I. ^ 0.556. 

^ The integrated intensity includes the three blended HPS components. 

Weakest HPS component at 93171.6210 M Hz w ith R.I. ^ 0.037 unless specified. At offsets (80",0'') we use the isolated HPS component at 
93176.2650 MHz with R.I. ^0.1111 (see Table fcTTl . 

^ A.V — W/T-i-nh is the equivalent width. Error bars are lo", and upper limits on W and 

/ Ntat X Nq/R.I. 

^ Fractional abundances with respect to H assuming N[H.) — 2N(}i2)- 
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Table D.2. As Table ID. II but for lines observed towards L 1544. 
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" Weakest HFS component at 113520.414 MHz with R.I. = 0.0184 

' For this offset, the strongest HFS component at 93173.7767 MHz with R.I.=0.2593 was used instead of the weakest. 
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Table D.3. As Table \DA\ but for lines observed towards Opli D. 
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Detections are at the level on W. 
^ Weakest HFS component at 88633.9360 MHz with R.I. ^ 0.111. 
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Table D.4. As Tabic ID. II but for lines observed towards L 1517B. 



Line 


6x 


5xj 


Tmb 


W 






^(H2) 


7Vtot/27V(H2) 




" 




mK 


mK km 


km ^ 


10^^ cm~^ 


10^2 cm~^ 






TTT" 

-10 


-20 


1 -, r-l 

115i 15 


r>ri_|_ 

30ih 3 


— — -— 

0.26 


— — 

0.88 


—3 

2.7 


—7 — —— 

1 .6( — 11) 




-50 


-20 


< 70 


< 20 




< 0.60 


1.8 


1 . 7( — 11) 




-ou 






24i 4 


n oi^ 


n 71 -1- n 1 o 


O T 


1 ^( 1 1 ^ 
1 .o^^ — i.i-) 




10 


-20 


1 nn-U on 


40di 5 


0.40 


1 .2ib 0.15 


1.9 


3. 1 ( — 11) 




30 


-20 


95i 25 


40ib 5 


0.42 


1 .2it 0.15 


1 . 1 


5.3( — 11) 




-10 


-60 


65ih 20 


30=h 5 


0.46 


n s Q -1- n 11^ 
U.ooit U.io 


1.5 


2.9( — 11) 




-10 


A n 


140ih 25 


28ih 3 


0.20 


n QO-I- n no 


2. 1 


2.0( — 11) 




1 n 

-lU 


n 
u 


1 9n-i- on 


OUi 4 


n oi^ 


n Hs-i- n 1 o 


2 2 


ri( 1 1 ^ 




-10 


20 


70± 20 


28± 5 


0.40 


f) QO+ n 1 


1.5 


2.7( — 11) 




-10 


-20 


280i 25 


llOih 5 


0.39 


0.38 


2.7 


7. 1 ( — 12) 




-ou 




1 fin-i- /in 




n "^1 


n 1 7-1- n n"? 


1 s 


A "?( 1 oA 

4. / ^^ — iz ) 




-30 


-20 


oTH-i- /in 


80=h 8 


0.35 


0.27 


2.3 


5.9( — 12) 




10 


-20 


oori-l- Q/l 


llOih 8 


0.38 


0.38 


1.9 


9.9( — 12) 




30 


-20 


265ih 37 


67ih 6 


0.25 


0.23 


1 . 1 


1 .o( — 11) 




-10 


-60 


onn-i- /I c: 


60ih 9 


0.30 


n on-i- n m 


1.5 


6.6( — 12) 




1 n 

-lU 




"^on-i- /in 

O-sUi 4U 


Q/1 -1- fi 

y4i o 


n OQ 
u.zy 


n '^o 


2 1 






-10 





295ih 35 


llOi 10 


0.37 


0.38 


2.2 


8. 7( — 12) 




-10 


20 


185± 40 


60± 8 


0.32 


0.20± 0.03 


1.5 


6.5( — 12) 




-10 


-20 


1105± 80 


683± 16 


0.62 


1.2 


2.7 


2.2(-ll) 




-50 


-20 


753± 71 


301± 12 


0.40 


0.52 


1.8 


1.4(-11) 




-30 


-20 


896± 76 


542 ± 15 


0.60 


0.94 


2.3 


2.0(-ll) 




10 


-20 


944± 71 


514± 14 


0.54 


0.89 


1.9 


2.3(-ll) 




30 


-20 


774± 73 


359± 13 


0.46 


0.62 


1.1 


2.7(-ll) 




-10 


-60 


661± 71 


390± 25 


0.59 


0.68 


1.5 


2.2(-ll) 




-10 


-40 


982± 75 


589± 15 


0.60 


1.0 


2.1 


2.4(-ll) 




-10 





1119± 70 


649± 14 


0.58 


1.1 


2.2 


2.5(-ll) 




-10 


20 


787± 71 


527± 15 


0.67 


0.91 


1.5 


2.9(-ll) 



Table D.5. As Table lD?3l but for lines observed towards L 310. The dust and thus H2 column densities are only the So- 
level. 



Line 


Sx 


Sy 


T,„b 


W 




JVtot 


N(H2) 


Aftot/2A'(H2) 






n 


niK 


niK km ^ 


km ^ 


lO^^cm-^ 


10^^ cm"^ 




CN 


30 


80 


390± 60 


230± 20 


0.5 


67 


1.5 


4.5(-9) 




10 


60 


< 305 


< 120 




< 35.2 


<0.7 






10 


80 


< 303 


< 130 




< 38.1 


<0.7 






30 


40 


< 288 


< 110 




< 32.2 


<0.7 






30 


60 


300± 55 


110± 15 


0.3 


32±4.4 


1.5 


2.1(-9) 




50 


60 


< 280 


< 110 




< 32.2 


<0.7 




HCN 


30 


80 


615± 35 


370± 10 


0.60 


6.1 


1.5 


4.1(-10) 




10 


60 


300± 55 


230± 20 


0.77 


3.8 


<0.7 






10 


80 


290± 60 


155± 20 


0.53 


2.5±0.33 


<0.7 






30 


40 


150± 50 


< 70 




< 1.15 


<0.7 






30 


60 


525± 35 


300± 15 


0.57 


4.9 


1.5 


3.3(-10) 




50 


60 


280± 55 


150± 25 


0.54 


2.5±0.41 


<0.7 






30 


80 


160± 20 


101± 7 


0.63 


0.34 


1.5 


2.3(-ll) 




-10 


80 


< 90 


< 36 




< 0.12 


<0.7 






10 


80 


< 90 


< 36 




< 0.12 


<0.7 






30 


40 


< 105 


< 42 




< 0.14 


<0.7 






30 


60 


150± 30 


70± 11 


0.47 


0.24±0.04 


1.5 


1.6(-11) 




30 


100 


< 95 


< 40 




< 0.14 


<0.7 






30 


120 


< 95 


< 40 




< 0.14 


<0.7 






50 


60 


< 135 


< 55 




< 0.19 


<0.7 






50 


80 


165± 35 


65± 8 


0.39 


0.22±0.03 


1.3 


1.7(-11) 




50 


100 


< 140 


< 60 




< 0.20 


1.5 


< 1.6(-11) 




70 


80 


< 90 


< 40 




< 0.14 


<0.7 






30 


80 


405± 25 


320± 10 


0.79 


0.55 


1.5 


3.7(-ll) 




-10 


80 


< 125 


< 55 




< 0.10 


<0.7 






10 


80 


205± 45 


100± 15 


0.49 


0.17±0.03 


<0.7 






30 


40 


< 120 


< 50 




< 0.09 


<0.7 






30 


60 


390± 40 


200± 10 


0.51 


0.35 


1.5 


2.3(-ll) 




30 


100 


150± 40 


70± 15 


0.47 


0.12±0.03 


<0.7 






30 


120 


< 120 


< 52 




< 0.09 


<0.7 






50 


60 


230± 60 


< 75 




< 0.13 


<0.7 






50 


80 


465± 40 


< 54 




< 0.09 


1.3 


< 6.9(-12) 




50 


100 


365± 55 


< 73 




< 0.13 


1.5 


< 8.7(-12) 




70 


80 


165± 40 


100± 10 


0.61 


0.17 


<0.7 





